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Abstract 

We perform a statistical study of permanent changes in longitudinal fields 
associated with solar flares by tracking magnetic features. The YAFTA feature 
tracking algorithm is applied to GONGH — h one-minute magnetograms for 77 
X- and M-class flares to analyze the evolution and interaction of the magnetic 
features and to estimate the amount of cancelled flux. We flnd significantly more 
magnetic flux decreases than increases occurred during the flares, consistent with 
collapsing loop structure. Correlations of both unsigned and signed flux changes 
with GOES peak X-ray flux are dominated by X-class flares in limb locations. 
The flux changes were accompanied in most cases by significant cancellation, 
most of which occurred during the flares. We find that the field strength and 
complexity at the polarity inversion line are approximately equally important 
in the fiux cancellation processes accompanied flares. We do not find correlation 
between the flux cancellation events and the stepwise changes in the features' 
fluxes. 
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1. Introduction 



It is known now that photospheric magnetic fields change significantly, abruptly 



and permanently as a result of major X- and M- class flares (e.g. Cameron 


and Sammisj 


1999| Kosovichev and Zharkova 2001 Wang et al.l 12002 |Sudol 


and Harveyl 


20051 Wang and Liu| |2010| |Petrie and Sudol, 2010; Wang et al. 


2012 


Gosain 


2012 Sun et al.| 2012| Petrie| 2012). Detection and study of pho- 



tospheric magnetic filed changes during solar flares can be important towards 
better understanding of the mechanisms of energy conservation and release in 
the flares. Despite the observations, most models of flares assume that these 
events are coronal and no photospheric magnetic fields changes occur across the 
flare duration (Priest and Forbes, 2002). Ideally, coronal field models should 
be based on boundary conditions observed in the lower atmosphere. However, 
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it is not clear how to relate approximately force-free current-free coronal field 
restructuring with the changes in the forced photosphere. 

A few possible mechanisms were offered in the past work that could cause 
sudden and permanent photospheric field changes observed. The most likely 
mentioned processes are the emergence of new magnetic flux and change in the 



direction of the field ( Spirock, Yurchyshyn and Wang, 2002 Sudol and Harvey, 



2005). Supporting magnetic flux emergence as a possible mechanism, Wang and 



Shi ( 1993 ) observed emergence of new flux in sunspot's umbra immediately after 
X-class flare. iFletcher and Hudson (2008) suggested the Alfven wave scenario as 
a possible explanation of the observed rapid variations to the LOS component 
of the photospheric magnetic field across the flare impulsive phase, when large- 
scale Alfven wave pulses transport energy and magnetic field changes rapidly 
from the flare site in the corona to the lower atmosphere. [Hudson, Fisher and| 



Welsch ( 2008 ) estimated back reaction of the photospheric magnetic field during 



a solar flare and concluded magnetic implosion and collapsing loop structure 
towards neutral lines. Wang and Liu (2010), [Petric and Sudol_(2010|, Wang et al. 



( 2012 ) found strong evidence that this prediction is true. Petrie and Sudol ( 2010 ) 
analyzed pixel-by-pixel the field changes accompanying 77 X- and M-class flares 
using GONG full-disk magnetograms. In this work, we use the flares from the 
GONG flare list to study flare-related photospheric field changes decomposing 
the field into individual magnetic features. 

As Sudol and Harvey^ {2005) have noted, the change in magnetic flux may 
be a more important physical quantity than the largest local field change at 
one particular location. Representative pixels generally represent the fastest and 



largest field changes, free of artifacts (Sudol and Harvey, 2005 Petrie and Sudol, 



2010), but do not represent a complete picture of the changes to the magnetic 



field in the region. On the other hand, magnetic flux calculations capture the 
effects of flares on entire active regions. The flux calculations of past authors 
have either summed the flux from changing pixels across entire active regions 
(Petrie and Sudol, 2010) or confined attention to simple bipolar active regions 



with strong neutral lines (e.g. Wang 



Wang 


2006 


Wang and Liu 


2010 



work the active regions are resolved into individual magnetic features whose flux 
changes are calculated individually. This allows us to identify which parts of a 
region change significantly and also allows us to investigate the interactions of 
different features in an active region. 

One of the aspects of the features' interaction associated with flares is a 
flux cancellation. Flux cancellation plays an important role in some theories of 
solar eruptions (e.g., van Ballegooijen and Martens 1989 Amari et al. 2010). 
It is found to be essential for the formation of prominences (Martin, 1998), and 



possible trigger of filaments eruptions and coronal mass ejections (Wang, Shi and 



Martin, 1996 Zhang et al., 2001 ). The flux cancellation mechanism is suggested 



to be a necessary condition of flare initiation as a part of slow reconnection 



processes in the lower atmosphere (Martin, Livi and Wang, 1985 Livi et al.. 



1989 Mathew and Ambastha, 2000 ). First evidence of the possibility of magnetic 



reconnection in flux cancellation was presented by Wang and Shi ( 1993 ) who 
analyzed vector magnetic field changes and demonstrated how features previ- 
ously not directly connected by transverse fields partially reconnected during 



SOLA: BurtsevaPetrie_2.tex; 17 October 2012; 0:24; p. 2 



flux cancellation. Strong observational evidence of magnetic reconnection in the 
lower atmosphere leading to instability and flare/CME events was presented by 



Zhang et al. 


(2001 


), 


Kim et al. 


(2001 


), 



Reconnection can occur at different levels of solar atmosphere: above, be- 



low or at the photospheric level (Zwaan, 1987). Harvey et al. (1999) observed 



magnetic flux disappeared in the chromosphere earlier than in the photosphere 
at the cancellation sites that suggests submergence of the flux implying that 
reconnection occurs above the photosphere. Chae (2004) confirms this finding. 



Yurchyshyn (2001), however, reported flux cancellation, followed a C-class flare, 



with the presence of upflows in the canceling feature, that suggest the reconnec- 
tion below or at level of the photosphere. In case of the prominences formation. 



cancellation is explained by emergence of ?7-shaped loops (Lites et al., 1995 Low, 



2001 Welsch, DeVore and Antiochos, 2005). Most of these studies suggest flux 



cancellation as pre-flare condition. Photospheric magnetic shear increases have 
been observed to occur during and after (Wang, 1992), and before and during 
(Wang et al., 1996) flares. In this work, we estimate flux cancellation in the 
features at the polarity-inversion lines, associated with the major flares, at the 
time before, during and after flare. 



2. Data analysis 



In this study full-disk photospheric magnetograms from GONG's sites are used. 
GONG full-disk magnetograms are obtained with 1-minute cadence, 2.5-arcsec 
pixel size, 3 G per pixel noise level. We have analyzed 77 sets of GONG magne- 
tograms for 38 X- and 39 M-class flares. We conflned our attention to the most 
energetic flares with the best data coverage. We eliminated all flares of GOES 
class weaker than M5.0 and all with an apparent central meridian longitude 
difference greater than 65 degrees. We further limited our attention to those 
flares for which GONG magnetograms are available of order one hour before and 
after the flare. The flares studied are identified in Tables 1 (M-class) and 2 (X- 
class) of Petrie and Sudol (2010). Petrie and Sudol ( 2010[ ) studied pixel-by-pixel 
the field changes accompanying the flares. In this work we study flare-related 
photospheric field changes analyzing collections of pixels (features). 



The preliminary data reduction procedure is as in Petrie and Sudol (2010). 



It is briefly described as follows: i) The images are remapped to azimuthal 
equidistant projection, 32° x 32° field of view, local to a flare, ii) The remaps 
are corrected for differential rotation and registered to a 10-minute average of 
magnetograms prior to the flare. 



We use the YAFTA feature-tracking code (Welsch and Longcope, 2003 ) to de- 



compose the field into individual magnetic features and to study their evolution 
and interplay during the flare. We have used YAFTA's "downhill" algorithm 
to group together more than 4 pixels lying on the same "hill" in absolute fleld 
strength above lOOG. The lOOG threshold identifles active region field in GONG 
magnetograms. Small features often have short lifetimes. The smallest feature 
that showed abrupt and permanent flux change in this work consisted of 8 pixels 
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Figure 1. Reference image, composed of 10 remapped GONG magnctograms, with tiie active 
region 10930 produced X6.5 fiare on 2006 December 6 and features id entified with tlie YAFTA 
downhill algorithm (left), c-map of stepwise field changes found by [Petrie and Sudol| | [2010[ l 
during this flare with locations of abrupt features' flux decreases (blue circles) and increases 
(red circles) found in this work over-plotted (right). 



(which is ~ 0.5 degrees in diameter), and the biggest one was 650 pixels (~ 3 de- 
grees in diameter). The program accomplishes matching of features between two 
consecutive images. Most of the features fragment, merge, appear, and disappear 
between steps. This makes the matching process complex. As a result, features 
change their labels. Only features that do not appear, disappear, fragment, or 
merge between steps "propagate" themselves, they keep their label through all 
steps. 

To track distinct features, for each flaring active region we identify the features 
on a reference image composed of the 10 images immediately preceding the flare. 
If some of these features-hills (two or more) appear as a singe feature- hill in most 
of the frames during the whole analyzed time interval and each of the individual 
features does not show significant temporal variations, they are merged into one 
feature on the reference image. Example of an active region with 18 features 
found is shown in Figure 1 (left). We locate each of these features on each of 
our images. We analyze four hours of data covering each flare. We calculate the 
total flux and the flux-weighted average x- and y-coordinates of each feature. 

Most previous feature-tracking work has applied a radial field correction to the 
magnetograms before processing, dividing the longitudinal fields by the cosine 
of the heliocentric angle (e.g. Welsch, Christe and McTiernan 2011). For most 



photospheric fields this is generally a valid assumption (Svalgaard, Duvall and 



Scherrer, 1978 Petrie and Patrikeeva, 2009 1 . However, for the photospheric fields 



in a flaring active region, particularly those fields associated with the strong, 
twisted magnetic structures directly involved in the eruption, this assumption 
may not be applicable. Petrie and Sudol (2010) reported evidence of this in the 
result that the average size of photospheric field change tends to increase as 
a function of distance from disk-center, indicating that the fields changing the 
most tend to be nearly horizontal. The radial field correction is necessary when 
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Figure 2. Evolution of the magnetic flux in one of the features recorded in AR10930 during 
an X6.5 flare on 2006 December 6. The three vertical lines indicate the GOES flare start, 
maximum and end times. 



features inside of an active region are tracked for long periods of time to ensure 



constant threshold during whole tracking interval (e.g., see Welsch, Christe and 



McTiernan 2011), while we track features only up to four hours. 



3. Flux changes 

We analyze size and sign of the flux change of each feature during the flares. The 
size and sign of the flux change are defined by fitting the total flux evolution in 
each feature before flare's onset and after flare end with straight lines assuming 
the flux evolution around flare is linear and ignoring long-term variations and 
trends. Figure 2 demonstrates an example of a strongest field change found in 
our analysis. The spread of the points is mainly caused by seeing changes in the 
Earth's atmosphere with some contribution from ongoing solar field changes. 

We selected features showing stepwise flux changes during flare. Each of the 
flux changes satisfies the following criterion. The value of the flux change is at 
least 1.3 a of the background flux scatter. The background flux is defined as a 
median value of the pre-flare feature's flux. The minimal and maximal values of 
the background fluxes are 2.4 x 10^^ and 7.3 x 10^^ Mx. The 1 a level is between 
1.8 X 10^^ and 4.5 x 10^° Mx. The maximal and median values of all the stepwise 
flux changes are 7.3 a and 3.0 a. It has been noted that duration of the flux 



changes is on the scale of 10-40 minutes (e.g. 


Kosovichev and Zharkova 


Wang and Liu 


2010 


Petrie and Sudol 2010 


). We select flux changes 



are complete within 30 minutes or less. The statistics for the flux increases 
and decreases for the 77 flares are presented in Table 1. In total 114 stepwise 
flux changes which satisfied the above criteria were selected. We find that field 
decreases are significantly more numerous than increases (70 versus 44, see first 
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Table 1. Number and absolute values of abrupt flux changes during flares 
(calculated by feature). The flux change values are in units of [10'^" Mx]. 



Field No. flux 
changes increases decreases 


Minimum 
change 


Maximum 
change 


Median 
change 


All 


44 


70 


0.02 


21.04 


0.61 


X-class 


21 


53 


0.02 


21.04 


0.81 


M-class 


23 


17 


0.10 


4.10 


0.48 


Disk-center 


15 


22 


0.10 


1.84 


0.35 


Limb 


29 


48 


0.02 


21.04 


0.82 


Disk-center X-class 


4 


11 


0.11 


1.84 


0.42 


Disk-center M-class 


11 


11 


0.10 


1.60 


0.27 


Limb X-class 


17 


42 


0.02 


21.04 


0.82 


Limb M-class 


12 


6 


0.15 


4.10 


0.84 



2 columns in line 1 of Table 1). The bias towards field decreases is noticeable for 
both disk-center and limb locations with slightly bigger ratio for the limb regions 
(first 2 columns in lines 4 and 5 of Table 1). We designate "near disk-center" for 
the locations within r < R/2 and "near limb" for the locations r > i?/2, where 
R is the solar disk radius in the image plane. Median longitudinal field changes 
are larger in the regions close to the limb (compare last 2 columns in lines 4 and 
5 of Table 1), indicating that their field vector is nearly horizontal in general. 
Limb changes are more than twice as numerous as disk center changes. X-flares 
show more numerous (74 versus 40) and strong photospheric flux changes in 
comparison to M-flares (lines 2 and 3 of Tabic 1). In case of X-class flares, flux 
decreases are more than twice as numerous as increases in both disk-center and 
limb regions, while for M-class flares ratio of decreases to increases is close to 
1:1 for disk-center and have the opposite bias 1:2 towards flux increases for close 
to limb regions. 

These statistics for magnetic features agree with the basic result reported by 



Petrie and Sudol (2010) that decreasing longitudinal fields tend to outnumber 



increasing longitudinal fields, except that here magnetic features of near-limb 
M-class flares have twice as many longitudinal flux increases as decreases. 

The overall bias towards the flux decreases in comparison with flux increases 
can be seen in the scatter plot of the change in the net magnetic flux against 
the background net magnetic flux shown Figure 3 (top left). The strongest flux 
change was found during XI. 2 flare on October 26, 2003 with onset at 05:57UT, 
which is much larger than all other flux changes found in this work. However, 
Figure 3 (top right), where the change in the unsigned magnetic flux against the 
background unsigned magnetic flux is plotted, demonstrates that this maximum 
flux change case is well aligned with the common flux change trend, showing 
that this value is not an outlier. This is a log-log plot. The point is of the same 
order of magnitude as the rest, but otherwise is an outlier. Most of the larger 
flux changes are associated with X-class flares. Figure 3 (bottom left), where the 
change in the unsigned magnetic flux against radial position on the solar disk is 
plotted, shows that the large X-class flares' flux changes occurred closer to the 
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Figure 3. Scatter plots of the change in the net magnetic flux against the background net 
magnetic flux (top- left), the change in the unsigned magnetic flux against the background 
unsigned magnetic flux (top-right) and against radial position on the solar disk (bottom-left) 
during M- (blue squares) and X-class (red circles) flares. In bottom-right panel, scatter plot of 
the change in the unsigned magnetic flux as a function of the GOES peak X-ray flux for active 
regions flared near disk center (red pluses) and near limb (blue crosses). The flux change is 
calculated by feature. The dotted lines at ±26 X 10^*^ Mx on the top-left plot separate features 
whose flux change negatively correlated with their background flux from the big features whose 
flux change did not conform the general pattern. 



limb. In total, a larger range of flux changes is found near the limb than near 
the disk center. The conclusions are consistent with the statistics of flux changes 
from [Petrie and Sudol (2010) (see Sections 5 and 7 of their paper). Figure 3 
(bottom right) , where the change in the unsigned magnetic flux as a function of 
the GOES peak X-ray flux is plotted, shows weak overall increasing trend of the 
flux change with flare strength. 

We note that features' sizes change causing an uncertainty in the flux esti- 
mates. Some fields crossing the 100 G threshold introduce this error into the 
flux estimate but such errors do not change the basic patterns of increasing and 
decreasing fluxes. 

We compute Pearson correlation coefficients between the magnetic fiux change 
and background magnetic flux. We find that flux change of the features with 
absolute value of background flux < 26 x 10^° Mx negatively correlated with 
their background flux, while the larger flux features did not conform the general 
pattern (see Figure 3, upper-left). The bimodal behavior seems to be analogous 
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Table 2. Pearson correlation coefficients for magnetic flux 
changes against background magnetic flux and GOES peak 
X-ray flux (calculated by feature). P- values are given in paren- 
thesis. 



Flux change vs. Flux change vs. 
background flux peak X-ray flux 



All 


-0.42 (6.6 X 10"'') 


0.38 (2.5 X IQ- 




X-class 


-0.50 (1.0 X 10-5) 


0.35 (2.4 X 10- 




M-class 


0.22 (0.19) 


0.07 (0.67) 




Disk-center 


-0.64 (1.5 X 10-5) 


0.13 (0.46) 




Limb 


-0.40 (4.1 X 10-"') 


0.36 (1.5 X 10" 




Disk-center X-class 


-0.74 (1.0 X 10-^) 


0.14 (0.64) 




Disk-center M-class 


0.10 (0.68) 


0.02 (0.93) 




Limb X-class 


-0.49 (1.3 X 10-"') 


0.33 (1.1 X 10- 




Limb M-class 


0.23 (0.37) 


-0.01 (0.96) 





to the different behavior of the weak/strong field changes < / > fOO G in 



Petrie and Sudol (20101. The correlation coefficients computed over features with 
absolute values of the background flux < 26 x 10^° Mx are shown in Table 2. 
P-value, the probability that the observed correlation occurred by chance, is also 
included in Table 2. We flnd modest negative correlation of the flux changes with 
the background flux for both disk-center and limb locations. The correlation is 
stronger for near disk-center locations and dominated by X-class flares. M-class 
flares show no correlation of the flux change against the background flux as in 
this case we find about equal number of the flux decreases and increases (see 
Table 1). 

Pearson correlation coefficients of the magnetic flux change against GOES 
peak X-ray flux show overall weak statistically significant correlation dominated 
by X-class flares in limb locations (see Table 2). In Table 2, the GOES peak 
X-ray flux and the position on the disk are approximately equally influential in 
the correlations. 

3.1. Flux changes calculated by flare 

The magnetic flux changes estimated in this work represent flux changes in a 
certain locations (features) . In order to have a fuller picture of flare effect on the 



active region, Petrie and Sudol (20101 summed the pixel-to-pixel field changes 
they found in each active region over the entire active region. We calculate 
flux change by flare summing features' abrupt flux changes during each flare in 
each active region. Some features in our analysis do not involve stepwise flux 
changes. Flux integrals over magnetic features may cause a mixture of some 
pixels that do not change value in an abrupt, stepwise way and some pixels with 
abrupt, stepwise changes of opposite sign. This can cause even features with 
many stepwise-changing pixels to have an overall flux profile that does not show 
a clear stepwise change. We restrict attention to those features whose overall flux 
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shows a clear, stepwise, permanent change. We expect these changes to dominate 
the flux evolution of an AR during a flare. 

An example of the c-map of stepwise fleld changes found by |Petrie and Sudol] 
(2010) during the X6.5 flare on December 6, 2006 with over-plotted locations of 
abrupt features' flux changes found in this work is shown in Figure 1 (right). Fea- 
ture 3, which demonstrated one of the strongest negative stepwise flux changes in 
this work (see Figure 2), is fully populated by pixels with strong stepwise changes 
of the same sign during the flare. Feature 1, which included both field decreases 
and increases and most of whose pixels did not exhibit clear, stepwise field 
changes because of spectral line changes during the flare, nevertheless showed a 
clear stepwise flux increase. Feature 5 also partially covered with field decreases, 
showed flux decrease, twice weaker than Feature 3. Feature 2 partially consisted 
of pixels with strongest negative field changes, but its fiux change did not exceed 
1.3 a of the background flux scatter criterion, as described above in this Section. 
The other features had smaller stepwise field decreases, not specified in this plot. 
They also showed stepwise flux changes twice smaller than Feature 3. For similar 
reasons, some of the 77 flares did not include a feature whose fiux changed in a 
clear stepwise manner even though each of these fiares did include at least one 
pixel that showed a stepwise change (Petrie and Sudol, 2010). In addition, the 
imposed field threshold restrict our analysis to strong field values only, which 
may also leave some weaker feature's flux changes beyond this analysis. 

Thus, from 77 flares analyzed only in 45 (25 X- and 20 M-class) fiares at 
least one feature with strong stepwise flux change was found. The active regions 
which did not involve stepwise flux changes during some flares do not seem to 
correlate with distance of the active region from the disk center. Some flares in 
the same active region and at similar separation from the disk center involved 
features with stepwise flux changes and others did not. Absence of the steps 
in the features fluxes also does not depend on a GONG observing station. One 
thing to note is that all cases of absence of the stepwise flux change are associated 
with the flare class XI. 8 and less, except one X3.8 flare on January 17, 2005, 
that occurred in disk-center location and did not involve stepwise flux changes. 
On the other hand, from 45 flares involved stepwise flux changes, only 12 were 
above flare class XI. 8. 

Net and unsigned flux changes statistics are shown in Tables 3 and 4, respec- 
tively. We find generally more decreases than increases in both net and unsigned 
flux changes (see the top line in Tables 3 and 4) . The bias is significant in case of 
the disk-center flares for both net and unsigned flux changes with corresponding 
ratios of the increases to decreases 2.6:1 and 2:1 (compare lines 4 in tables 3 and 
4) . Limb flares show approximately equal number of increases and decreases for 
both net and unsigned flux changes (lines 5 in tables 3 and 4). 



Petrie and Sudol (2010) plotted the longitudinal field components associated 



with possible magnetic vector tilt and azimuth angles at different locations on the 
solar disk. They found that if field vectors become more horizontal during every 
fiare then the longitudinal fields can either increase or decrease, both near disk 
center or near the limb. However, decreases would likely outnumber increases at 
all parts of the disk investigated (i.e., within 65 degrees of disk center) and this 
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Table 3. Net flux changes during flares (calculated by flare). The flux change 
values are in units of [10^'' Mx]. 



Field No. fleld 
changes increases decreases 


Minimum 
change 


Maximum 
change 


Median 
change 


All 


18 


27 


0.002 


22.17 


1.17 


X-class 


8 


17 


0.002 


22.17 


1.63 


M-class 


10 


10 


0.12 


2.52 


0.80 


Disk-center 


5 


13 


0.002 


2.65 


0.63 


Limb 


13 


14 


0.29 


22.17 


1.75 


Disk-center X-class 


2 


5 


0.002 


2.65 


1.32 


Disk-center M-class 


3 


8 


0.12 


2.07 


0.55 


Limb X-class 


6 


12 


0.29 


22.17 


2.89 


Limb M-class 


7 


2 


0.30 


2.52 


1.36 



Table 4. Unsigned flux changes during flares (calculated by flare). The flux 
change values are in units of [10^'^ Mx]. 



Field No. fleld 
changes increases decreases 


Minimum 
change 


Maximum 
change 


Median 
change 


All 


19 


26 


0.04 


17.50 


1.16 


X-class 


6 


19 


0.29 


17.50 


1.75 


M-class 


13 


7 


0.04 


2.16 


0.63 


Disk-center 


6 


12 


0.07 


2.65 


0.84 


Limb 


13 


14 


0.04 


17.50 


1.75 


Disk-center X-class 


1 


6 


0.37 


2.65 


1.32 


Disk-center M-class 


5 


6 


0.07 


2.16 


0.55 


Limb X-class 


5 


13 


0.29 


17.50 


2.73 


Limb M-class 


8 


1 


0.04 


2.14 


1.36 



bias is expected to be greater near disk center than near the hmb. We find these 
patterns in our data. 

For the X-class flares, decreases are at least twice as numerous as increases in 
both disk center and limb locations and for both net and unsigned flux changes 
(see lines 2, 6 and 8 in tables 3 and 4). Similar bias is reported by Petrie and 



Sudol (20101 for the X-class unsigned flux changes while for the net flux changes 



they do not find any significant bias. We find more than twice more decreases 
than increases for the M-class net flux changes in disk-center locations, which 
is consistent with the same statistics from Petrie and Sudol (20101, while the 
bias switches to the opposite for the limb flares (compare lines 3, 7 and 9 in 
table 3). For the M-class unsigned flux changes, the increases are nearly twice as 
numerous as decreases in general (lines 3 in table 4). A significant bias towards 
increases in near-limb M-class unsigned flux is observed, while same category 
from Petrie and Sudol (20101 does not show any significant pattern. Similar to 
the feature-wise statistics, median flux changes are larger in the regions close to 
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Figure 4. Scatter plots of the change in the net magnetic flux and unsigned magnetic flux 
against the background net and unsigned magnetic flux, respectively, (top row), against radial 
position on the solar disk (middle row) during M- (blue squares) and X-class (red circles) flares, 
and as a function of the GOES peak X-ray flux for active regions flared near disk center (red 
pluses) and near limb (blue crosses) in the bottom-row plots. The flux change is calculated by 
flare. 



the limb for both unsigned and net flux changes and larger in case of X-flares in 
comparison with M-class. 

Figure 4 shows scatter plots of the change in the net and unsigned magnetic 
flux against the background net and unsigned magnetic flux, against radial 
position on the solar disk as well as a function of the GOES peak X-ray flux. 
The correlations of magnetic flux changes against background magnetic flux and 
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Table 5. Pearson correlation coefficients for magnetic flux 
changes against background magnetic flux (calculated by flare) . 
P-values are given in parenthesis. 



Background flux vs. 





net flux 


unsigned flux 




change 


change 


All 


-0.28 (0.09) 


0.55 (7.0 X 10-5) 


X-class 


-0.39 (0.08) 


0.55 (4.0 X 10-3) 


M-class 


0.32 (0.20) 


0.47 (3.4 X 10-2) 


Disk-center 


-0.56 (2.4 X 10-2) 


0.68 (1.0 X 10-3) 


Limb 


-0.26 (0.23) 


0.56 (2.3 X 10-3) 


Disk-center X-class 


-0.63 (0.19) 


0.68 (0.10) 


Disk-center M-class 


-0.42 (0.24) 


0.87 (5.5 X 10-5) 


Limb X-class 


-0.41 (0.13) 


0.52 (2.5 X 10-2) 


Limb M-class 


0.55 (0.17) 


-0.33 (0.44) 



Table 6. Pearson correlation coefficients for magnetic flux 
changes against GOES peak X-ray flux (calculated by flare). 
P-values are given in parenthesis. 



Peak X-ray flux vs. 





net flux 


unsigned flux 




change 


change 


All 


0.49 (6.5 X 10-"') 


0.56 (8.5 X 10-5) 


X-class 


0.39 (0.07) 


0.47 (1.8 X 10-2) 


M-class 


-0.11 (0.65) 


0.21 (0.40) 


Disk-center 


0.34 (0.17) 


0.36 (0.14) 


Limb 


0.48 (1.4 X 10-2) 


0.56 (3.9 X 10-3) 


Disk-center X-class 


0.20 (0.72) 


0.20 (0.68) 


Disk-center M-class 


-0.15 (0.65) 


0.09 (0.80) 


Limb X-class 


0.40 (0.11) 


0.50 (4.1 X 10-2) 


Limb M-class 


0.003 (0.99) 


0.81 (2.2 X 10-2) 



GOES peak X-ray flux are summarized in Tables 5 and 6. The only statistically 
significant modest negative correlation of the net flux changes against net back- 
ground flux found is for the disk-center flares (see first column in line 4 of Table 
5), which agrees with the results in Table 3. Some weak negative correlation 
can also be noted in case of the X-class flares (first column in line 2 of Table 
5), which is consistent with the top- left Figure 4 and statistics in Table 3. The 
distribution of points in the figure is not much different from the one shown 
in Figure 3. The top-right Figure 4 shows more points that do not conform 
the general trend for the M-class unsigned fiux changes, while the X-class flux 
changes maintain the pattern quite well, except one point, corresponding to one 
of the limb flares. The middle-row plots in Figure 4 clearly show that stronger 
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flux changes occur preferentially closer to the limb. Also a larger range of the 
flux changes corresponds to the limb locations. 

The bottom-row plots in Figure 4 show overall increasing trend of the net and 
unsigned flux changes as a function of the peak X-ray flux. We do find overall 
moderate statistically significant correlation with the peak X-ray flux for both 
net and unsigned flux changes dominated by the X-class flares in limb locations, 
that is more clear in case of unsigned flux change, while the net flux change 
shows statistically significant contribution to the linear correlation only for limb 
flares (compare lines 1, 2, 5 and 8 of Table 6). 



4. Flux cancellation 



We also identify episodes of flux cancellation, when closely-spaced opposite- 
polarity features simultaneously mutually lose flux (Livi, Wang and Martin, 



1985 Martin, Livi and Wang, 1985). We estimate the flux cancellation as de- 



scribed in Welsch, Christe and McTiernan (2011 ). Features' flux computed from 



high-cadence data suffer short-term fluctuations (noise) with large amplitudes 
in the case of this study as GONG data also influenced by earth's atmosphere. 
The fluctuations significantly affect flux cancellation measurements leading to 
exaggerated values of the flux losses in some cases while summing up noise 
instead of real flux evolution. We smooth the total flux in each feature with 30- 
minute running mean. After that we mostly have long-term evolutionary trend 
of the flux that we use in further flux cancellation computation. First, we de- 
fine instances when boundaries of closely-spaced opposite-polarity features come 



into contact, using the technique described in [Welsch, Christe and McTiernan 
(2011). In case of our analysis, the closest distance between features, confined 
by GONG's pixel-size, is about 1.8 Mm on the Sun at disk center. Further, for 
a pair of neighboring features to be partially cancelled over a time step, the 
features centers of flux are required to approach and the fiux of both features to 
decrease. This approach is subject to many sources of uncertainties, related to 
threshold applied (some flux is not included) , fiuxes fluctuations from noise and 
true evolution, as well as possible interaction with several neighboring like- or 
opposite-polarity features. To account for them, flux cancellation is computed 



in three ways as in Welsch, Christe and McTiernan (2011 ): 



minimum flux loss from the pair of canceling features, summed over all 
cancellation events; 

average flux loss from the pair of features, observed to cancel over at least 
two consecutive steps (multi-step cancellation), summed over all multi-step 
cancellation events; 

minimum flux loss from the pair of features, observed to cancel over at least 
two consecutive steps, summed over all multi-step cancellation events. 

The total number of cancellation events and the estimates of the flux losses 
during the whole 4-hour interval around each flare for M- and X-class flares are 
listed in Tables 7 and 8, respectively. Ordering of all the estimates obtained with 



the three methods for each flare is consistent with that from Welsch, Christe and 
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Table 7. Cancelled flux values (in units of [10^'' Mx]) for 4- hour intervals around each M-class flare, computedl 
in three ways (see text). | 



Date GOES Radial position, No. No. multi- All events. Multi-step Multi-step 





class 


in Rq 


events 


step events 


J2 min{5^) 


J2avg{5<!>) 


inin{5^) 


2001 Jun 22 


M6.2 


0.74 












2001 Jun 23 


M5.6 


0.43 


29 


2 


5.5 


0.2 


0.2 


2001 Sep 05 


M6.0 


0.49 


9 




0.7 






2001 Sep 09 


M9.5 


0.70 


14 




1.1 






2001 Sep 16 


M5.6 


0.90 


2 




0.4 






2001 Oct 22 


M6.7 


0.52 


35 


8 


2.2 


1.3 


0.6 


2001 Oct 23 


M6.5 


0.43 


39 


13 


5.0 


5.7 


2.5 


2001 Nov 07 


M5.7 


0.74 


96 


21 


18.2 


7.1 


4.6 


2001 Nov 08 


M9.1 


0.48 


5 




0.4 






2001 Nov 28 


M6.9 


0.28 


53 


13 


7.7 


3.1 


2.0 


2001 Nov 29 


M5.5 


0.19 


133 


29 


11.2 


4.2 


2.5 


2001 Dec 26 


M7.1 


0.83 


26 


4 


12.9 


4.8 


3.6 


2002 Jan 09 


M9.5 


0.29 


60 


9 


10.8 


3.7 


2.4 


2002 Mar 14 


M5.7 


0.41 


28 


9 


3.6 


2.7 


1.5 


2002 Jul 11 


M5.8 


0.80 












2002 Jul 17 


M8.5 


0.33 


26 


7 


1.1 


1.7 


0.5 


2002 Jul 26 


M8.7 


0.49 


173 


38 


21.7 


9.5 


6.2 


2002 Aug 16 


M5.2 


0.39 


62 


11 


11.7 


4.4 


3.0 


2002 Oct 05 


M5.9 


0.52 


10 


2 


0.3 


0.2 


0.1 


2002 Nov 18 


M7.4 


0.80 












2002 Dec 20 


M6.8 


0.65 


183 


52 


12.2 


7.1 


3.9 


2003 Oct 26 


M7.6 


0.63 


95 


12 


44.8 


13.3 


6.1 


2003 Oct 27 


M5.0 


0.54 


110 


28 


11.0 


10.8 


3.3 


2003 Nov 20 


M5.8 


0.33 


45 


9 


3.5 


1.5 


0.5 


2003 Nov 20 


M9.6 


0.18 


88 


28 


7.9 


6.2 


2.6 


2004 Jan 17 


M5.0 


0.37 












2004 Jan 20 


M6.1 


0.32 


7 




0.1 






2004 Jul 13 


M6.2 


0.83 


10 




0.8 






2004 Jul 13 


M6.7 


0.72 


14 


2 


2.5 


0.5 


0.5 


2004 Jul 20 


M8.6 


0.58 












2004 Jul 22 


M9.1 


0.35 


31 


2 


1.3 


0.04 


0.01 


2004 Jul 25 


M7.1 


0.39 


21 


6 


1.5 


1.1 


0.4 


2004 Aug 14 


M7.4 


0.56 


43 


10 


8.5 


2.5 


1.3 


2004 Oct 30 


M5.9 


0.39 


37 


10 


3.1 


4.9 


0.6 


2005 Jan 15 


M8.4 


0.33 


62 


12 


16.9 


5.7 


2.0 


2005 Jan 15 


M8.6 


0.31 


63 


12 


20.8 


10.7 


5.3 


2006 Dec 06 


M6.0 


0.90 












2007 Jun 04 


M8.9 


0.75 


34 


7 


2.4 


1.5 


0.5 
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Table 8. Cancelled flux values (in units of [10^^ Mx]) for 4-hour intervals around each X-class flare,! 
computed in three ways (see text). | 



Date GOES 
class 


Radial position, 
in Rq 


No. 

events 


No. multi- 
step events 


All events, 
J2 min{5^) 


Multi-step 


Multi-step 
J] mm(5<I>) 


2001 Apr 02 X20.0 


0.93 


54 


17 


63.5 


23.9 


17.0 


2001 Jun 23 XI. 2 


0.41 


42 


12 


5.9 


2.3 


1.3 


2001 Aug 25 X5.3 


0.66 


58 


10 


2.7 


0.9 


0.4 


2001 Oct 19 XI. 6 


0.50 


27 


7 


1.8 


0.8 


0.5 


2001 Oct 22 XI. 2 


0.47 


26 


9 


5.8 


4.8 


2.4 


2001 Dec 11 X2.8 


0.65 


8 


2 


0.8 


0.9 


0.3 


2002 May 20 X2.1 


0.91 


2 




0.2 


" 


" 


2002 Aug 21 XI. 


0.72 


19 


■ 


5.9 


3.8 


3.0 


2003 Mar 17 XI. 5 


0.59 


70 


26 


28.6 


17.6 


12.5 


2003 Mar 18 XI. 5 


0.71 


113 


31 


24.2 


10.7 


6.6 


2003 May 27 XI. 3 


0.36 


72 


11 


11.2 


3.8 


2.4 


2003 May 28 X3.6 


0.38 


84 


17 


12.1 


6.1 


3.0 


2003 Jun 10 XI. 3 


0.66 


96 


14 


11.7 


3.7 


1.7 


2003 Jun 11 XI. 6 


0.85 


77 


11 


14.5 


3.2 


1.6 


2003 Oct 19 XI. 1 


0.81 


14 




4.5 


" 




2003 Oct 26 XI. 2a 


0.73 


84 




25.1 


16.4 


■ 


2003 Oct 26 XI. 2b 


0.62 


127 


■ 


33.3 


15.2 


4.0 


2003 Oct 29 XIO.O 


0.34 


153 


37 


40.4 


25.6 


14.4 


2003 Nov 02 X8.3 


0.85 


" 




" 


" 


" 


2004 Feb 26 XI. 1 


0.43 


19 




2.7 


0.8 


0.4 


2004 Jul 15 XI. 6 


0.72 


44 


■ 


12.4 


2.8 


1.5 


2004 Jul 15 XI. 8 


0.79 


2 




0.2 


" 


" 


2004 Jul 16 XI. 1 


0.56 


59 




14.0 


2.4 


1.4 


2004 Jul 16 XI. 3 


0.63 


51 


■ 


10.6 


6.4 


3.2 


2004 Jul 16 X3.6 


0.56 


71 


16 


18.7 


6.0 


3.7 


2004 Aug 13 XI. 


0.49 


63 


16 


4.9 


1.9 


1.4 


2004 Oct 30 XI. 2 


0.40 


29 


4 


2.2 


0.4 


0.1 


2005 Jan 01 XI. 7 


0.56 












2005 Jan 15 XI. 2 


0.32 


194 


35 


16.0 


7.5 


2.7 


2005 Jan 15 X2.6 


0.30 


297 


95 


28.3 


19.2 


7.8 


2005 Jan 17 X3.8 


0.46 


203 


54 


73.2 


31.2 


17.4 


2005 Jan 20 X7.1 


0.88 


7 




7.6 






2005 Jul 30 XI. 3 


0.86 












2005 Sep 10 XI. 1 


0.75 


39 


3 


9.8 


3.2 


0.6 


2005 Sep 10 X2.1 


0.73 


43 


8 


12.5 


4.7 


1.9 


2005 Sep 13 XI. 5 


0.32 


120 


15 


9.5 


2.2 


1.1 


2005 Sep 15 XI. 1 


0.40 


90 


17 


16.3 


5.6 


3.7 


2006 Dec 06 X6.5 


0.89 


60 


10 


44.1 


14.4 


8.0 


2006 Dec 14 XI. 5 


0.75 


48 


8 


7.5 


4.5 


1.8 
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Figure 5. Estimates of flux [in Mx] canceled in the active regions over 60 minutes before 
(blue), 60 minutes after (red) and during (green) flare in multi-step events when flux can- 
cellation observed over at least two consecutive steps. Time during flare varies from 8 to 57 
minutes. Two top a- and two bottom b-panels show results for M- and X-flares, respectively. 
In case of each flare class: top panel shows sum of minimum flux loss for each pair of features 
in each event; bottom panel shows the total flux loss normalized by time interval. Flux losses 
are plotted in ascending order of flare class. The inset plots show the total flux loss over the 
entire time interval around each flare (one bar per flare). 
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I I 

McTiernan (2011 ). The third technique is the most restrictive, thus it should be 



the most reUable, showing real physical pattern. 

We also compute flux cancellations separately for three time intervals: before 
flare's onset, during flare and after flare end. The GOES flare start, pick and 
end times are used. By straight division of our data sets by three time intervals 
we end up with different duration of each of the intervals for each flare. For fair 
comparison, we restrict time intervals before and after flare to 60 minutes each. 
In case if a flare's data set is shorter at least in one of the two time intervals, 
the flare is excluded from the comparison. 29 X- and 29 M-class flares are flnally 
included. Time during flare is different for all flares and varies from 8 to 57 
minutes. Note, that the during-flare interval is always shorter than 60 minutes. 
Flux losses for 29 X- and 29 M-class flares obtained by the third, the most 
restrictive, technique for estimating flux cancellation are shown in Figure 5. We 
plot both the total flux losses and the flux losses normalized by time interval. 
We compute mean and median values of flux losses in flux cancellation events 
before, during and after the flare for both the total flux loss and the flux loss 
per minute. The results are shown in Table 9. 

We flnd that the amount of flux cancelled was greater over the time during the 
flare than before and after the flare for X-class flares. The ratios of the flux losses 
before, during and after the flare to the amount of flux cancelled before the flare 
are 1:3:0.7 (for median values) and 1:4:0.5 (mean values). Some of the strongest 
X-class flares showed significantly greater flux loss during the flare time. Flux 
loss before the flare was greater than after the flare by a smaller factor. In the 
case of the M-class flares no significant bias was found. 

The flux losses in Figure 5 are plotted in the ascending order of flare class. 
We note an overall increasing trend of the flux losses as a function of flare power. 
The trend is more obvious for X-class flares than for M-class flares. It is well 
seen in the inset plots of Figure 5, where each bar represents flux loss over the 
entire 4-hour interval around each flare (from last column of Tables 7 and 8) also 
plotted in the ascending order of flare class. Can the correlation of the flux loss 
with total flux in the canceling features or total flux in active region account for 
this trend? 



Table 9. Statistics of flux losses in flux cancellation events before, during and after 
flare for each of the computation techniques shown in Figure 5. 







Before flare 


During flare 


After flare 




Total flux loss in M-flares 


5.2 


1.8 


6.6 


Median flux loss 


X-flares 


6.1 


7.1 


4.4 


[1018 Mx] 


Flux / minute in M-flares 


0.09 


0.1 


0.1 




X-flares 


0.1 


0.3 


0.07 




Total flux loss in M-flares 


5.5 


6.5 


9.4 


Mean flux loss 


X-flares 


12.9 


22.3 


6.3 


[IQlS Mx] 


Flux / minute in M-flarcs 


0.09 


0.2 


0.16 




X-flarcs 


0.2 


0.8 


0.1 
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4.1. Correlation of the flux loss with X-ray flux. 



Flux loss, total flux in active region, number of canceling features and total flux in 
the canceling features against GOES peak X-ray flux are shown in Figure 6. The 
Pearson correlation coefficients for all the quantities are summarized in Table 10. 
The increasing trend of the flux loss as a function of flare strength dominated by 
X-class flares can be seen in Figure 6 for the flux loss and the total flux in active 
region, but less obvious for the number of canceling features and the total flux 
in canceling features. We find that all four quantities show weak-to-moderate 
positive correlation with flare strength in disk-center locations dominated by 
X-class flares, except the number of canceling features which shows only overall 
correlation with the peak X-ray flux for disk-center flares, but no statistically 
significant correlation for the disk-center X-class flares (lines 4 and 6 of Table 10). 
The correlation is highest for the flux loss, which also shows moderate positive 
correlation for the limb flares dominated by the X-class flares in limb locations 
(first column of Table 10). The M-class flares reveled no correlation with the 
peak X-ray flux in any of the categories. 
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Figure 6. Flux loss, total flux in active region, number of canceling features and total flux in 
the canceling features against GOES peak X-ray flux for active regions flared near disk center 
(red pluses) and near limb (blue crosses). 



In previous studies, it has been found that the total flux in active region and 



flux at the polarity-inversion line (PIL) correlate with flare activity (Leka and 
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Table 10. Pearson correlation coefficients for flux loss, total flux in active region, number of cancelingl 
features and total flux in the canceling features against GOES peak X-ray flux for active regions flared 
near disk center and near limb. P-values are given in parenthesis. | 
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active region 


features 


canceling features 


AU 


0.66 (5.8 X 10- 




0.07 (0.57) 


0.08 (0.56) 


0.09 (0.49) 


X-class 


0.66 (1.8 X 10- 




0.16 (0.32) 


0.01 (0.95) 


0.10 (0.60) 


M-class 


0.12 (0.58) 




0.18 (0.29) 


-0.02 (0.93) 


-0.02 (0.93) 


Disk-center 


0.74 (1.3 X 10- 




0.35 (4.8 X 10-2) 


0.47 (8.7 X 10-3) 


0.56 (1.2 X 10-3) 


Limb 


0.68 (2.8 X 10- 




-0.01 (0.94) 


-0.19 (0.33) 


-0.11 (0.56) 


Disk-center X-class 


0.75 (2.2 X 10- 




0.61 (2.5 X 10-2) 


0.39 (0.20) 


0.64 (1.6 X 10-2) 


Disk-center M-class 


0.28 (0.30) 




0.27 (0.25) 


0.24 (0.37) 


0.22 (0.43) 


Limb X-class 


0.70 (5.5 X 10- 




0.05 (0.81) 


-0.22 (0.38) 


-0.12 (0.63) 


Limb M-class 


-0.14 (0.71) 




0.02 (0.92) 


-0.44 (0.21) 


-0.32 (0.38) 



Barnes, 2007 Schrijver, 2007 Barnes and Leka, 2008 Welsch et al., 20091. The 



latter was found to correlate with peak X-ray flux the most (Schrijver, 2007). We 
find about the same level of low-to-moderate correlation of both the flux in the 
canceling features and the total flux in active region with flare strength for the 
disk-center locations dominated by X-class flares. Here, the canceling features 
are only the features contributed to the multi-step flux cancellation events, thus 
the flux in the canceling features and the flux at the PIL are not exactly the 
same quantites. The relation of the flux in PIL features to flare strength will be 
discussed in the next Section. 

We have explored how the flux loss, total flux in active region, number of 
canceling features and total flux in the canceling features correlate with each 
other. The results are shown in Table 11. We find that both the number of 
canceling features and the total fiux in the canceling features show moderate 
positive correlation with the fiux loss in all categories except limb flares, and in 
particular limb X-class flares, which did not show any correlation with the flux 
loss (last two columns in Table 11). However, some moderate correlation is found 
for the limb M-class flares (see line 9 of the last two columns in Table 11). The 
correlation for the flux in the canceling features is dominated by the disk-center 
X-class flares. In contrast, the total flux in active region show weak-to-moderate 
correlation only with the flux loss for the disk-center locations and with only 
contribution of the X-class flares (first column in Table 11). We also note that 
the flux loss is approximately equally correlated with the number of canceling 
features and the total flux in these features (see last two columns in Table 11). 
This points out that the field strength and complexity are approximately equally 
important. 

Thus, the correlation level of flux loss with the total flux in the canceling 
features is approximately as high as the correlation of the flux loss with the flare 
strength and can account for the correlation of the flux loss with the GOES 
peak X-ray flux we observe, in general, and in particular for the X-class flares in 
disk-center locations. However, we do not flnd correlation of the flux loss with 
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the total flux in the cancehng features (neither with the number of canceling 
features) for the limb X-class flares, thus it does not explain the correlation of 
the flux loss with the peak X-ray flux in this case. The correlation appears to 
work better at disk-center locations than far from disk center suggesting that 
it is the radial flux loss that correlates best with peak X-ray flux. We need 
vector data to know the radial flux and not the nearly horizontal field seen in 
line-of-sight images near the limb to confirm this conjecture. 



4.2. Correlation of the fiux at the PIL with X-ray flux. 



Schrijver ( 2007 ) measured unsigned flux near the polarity-inversion lines, R, and 



proposed this quantity as a measure for major flare potential rather than the 
total flux in an active region. In this work, the total flux in PIL features can be 
considered as an equivalent of the parameter R. Figure 7 shows total flux in the 
PIL features against total flux in the active region, as well as the scatter plot 
of the GOES peak X-ray flux against total flux in the PIL features. Total flux 
in the PIL features correlated with total flux in the active region with overall 
correlation coefficient ^ 0.60, which is in agreement with Schrijver (2007). The 
right plot indicates that large flares associated with strong PIL for the disk-center 
locations, consistent with the finding of I Schrijver (2007). 

Such strong-gradient PILs are generally associated with flux emergence and 
flux cancellation driven by photospheric flows. Total evolutionary net flux change 
in the PIL features during time before and after flare shows that PIL flux 
decreases before flare outnumber increases in X-class flares, while M-class flares 
show approximately equal number of flares with the increases and decreases. This 
could hint that there is more flux cancelled than emerged on PILs before flare. 
During time after flare we find approximately equal number of flux decreases 
and increases in both X- and M-class flare PILs. 



Table 11. Pearson correlation coefficients for flux loss against total flux in active 
region, against number of canceling features and against total flux in the canceling 
features for M- and X-class flares occurred near disk center and near limb. P-values 
are given in parenthesis. 



flux loss 


total flux in 


No. canceling 


r 


total flux in 






active region 


features 




canceling features 


AU 


0.21 (0.11) 


0.47 (1.4 X IQ- 




0.50 (5.5 X 10- 




X-class 


0.29 (0.10) 


0.40 (2.1 X 10" 




0.51 (2.2 X 10- 




M-class 


0.31 (0.12) 


0.70 (3.8 X 10- 




0.60 (9.8 X 10- 




Disk-center 


0.40 (3.3 X 10-2) 


0.70 (1.2 X 10- 




0.78 (8.8 X 10- 




Limb 


0.004 (0.98) 


0.14 (0.47) 




0.21 (0.27) 




Disk-center X-class 


0.63 (1.9 X 10-2) 


0.66 (1.2 X 10- 




0.86 (4.4 X 10- 




Disk-center M-class 


0.25 (0.36) 


0.69 (2.2 X 10- 




0.58 (1.6 X 10- 




Limb X-class 


-0.10 (0.69) 


0.02 (0.93) 




0.15 (0.56) 




Limb M-class 


0.41 (0.25) 


0.69 (2.4 X 10- 


^) 


0.61 (0.06) 
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Figure 7. Left: Total flux in the PIL features against total flux in active regions flared near 
disk center (red pluses) and near limb (blue crosses). Right: Total flux in the PIL features 
against GOES peak X-ray flux. The peak X-ray flux of 1 corresponds to XI class flare. 

4.3. Correlation of the flux loss with radial position on the solar disk. 

Figure 8 shows flux loss, total flux in active region, number of canceling features 
and total flux in the canceling features against radial position on the solar disk. 
The Figure suggests that the flux loss, total flux in active region and total flux in 
canceling features are not correlated with radial position on the disk. These fields 
have no apparent directional preference, whereas middle-row plots in Figure 4 
show correlation of the flux with radial position consistent with nearly horizontal 
fields. On the other hand, the bottom-left plot of Figure 8 does show correlation, 
indicating that the number of canceling features per active region is affected by 
foreshortening, decreasing towards the limb. Vector data is needed to clarify 
whether the heliographic spatial resolution is the reason for this bias. 

The Pearson correlation coefficients, summarized in Table 12, show overall 
weak negative correlation of the total flux in active region with radial position. 
The number of canceling features indicate slightly higher, weak-to-moderate 
negative correlation with radial position for the X-class flares only. The flux 
loss and flux in the canceling features show no statistically significant negative 
correlation with radial position. 



Table 12. Pearson correlation coefflcients for flux loss, total flux in active region, 
number of canceling features and total flux in the canceling features against radial 
position on the solar disk for M- and X-class flares. P-values are given in parenthesis. 





flux loss 


total flux in 


No. canceling 


total flux in 






active region 


features 


canceling features 


All 


0.18 (0.19) 


-0.23 (4.0 X 10~2) 


-0.18 (0.19) 


-0.08 (0.53) 


X-class 


0.10 (0.61) 


-0.24 (0.14) 


-0.43 (4.8 X 10^2) 


-0.18 (0.33) 


M-class 


0.17 (0.40) 


-0.21 (0.22) 


0.05 (0.80) 


-0.02 (0.93) 
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Figure 8. Flux loss, total flux in active region, number of canceling features and total flux in 
the canceling features against radial position on the solar disk. 



Cancellation events were not observed in all 77 flares . In general, absence of 
flux cancellation events correspond to the active regions flared far from the disk 
center where longitudinal component of magnetic field is decreased in comparison 
to disk-center locations, while radial component of magnetic field most likely 
would involve some cancellation events. However, there are exclusions when disk 
center flares do not involve flux cancellation or near limb flares show significant 
flux loss. 

The M5.0 flare on January 17, 2004, occurred at ~ 0.37 solar radius in 
NOAA10540, did not involve flux cancellation events. The active region is diffuse 
with total flux in all features lower than median value of total fluxes of the all 
active regions analyzed. The same active region hosted another flare from our 
analysis, M6.1 on January 20, 2004, which also occurred at near disk-center lo- 
cation (~ 0.32) and showed very small amount of flux loss with no multistep flux 
cancellations. The active region location coincides with location of NOAA10486, 
associated with Halloween's flares, only three solar rotations earlier. After two 
following rotations (NOAA10508 and 10523) the active region almost completely 
decayed. In another rotation, it appeared on the front side with new flux emerged 
into it as NOAA10540, hosting a few C- and M-class flares. 

Three other flares at near disk-center locations (~ 0.48, 0.35 and 0.49) which 
did not involve multistep flux cancellations (M9.1 flare on November 8, 2001, 
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M9.1 flare on July 22, 2004 and M6.0 flare on September 5, 2001), had diffuse 
structure of the active region with short PIL, in case of the first flare, while in 
the case of the last two flares, concentrated large total flux active region with 
total flux in the PIL features lower than median value of that for the all active 
regions analyzed. 

All active regions hosted X-class flares in the near disk-center locations in- 
volved flux cancellation events. Absence of flux cancellation events around X- 
class flares generally coincided with far from disk-center location of the flared 
active region. We note the opposite exclusions, when some active regions, flared 
in close to limb locations, involved significant flux cancellations during M- and 
X-class fiares (see Tables 7 and 8 for details): X20.0 fiarc on April 2, 2001 
(NOAA9393), M5.7 on November 7, 2001 (NOAA9690); M7.1 on December 
26, 2001 (NOAA9742); M6.8 on December 20, 2002 (NOAA10226); M7.6 on 
October 26, 2003 (NOAA10484); XI. 5 on March 17, 2003 (NOAA10314): X6.5 
on December 6, 2006 (NOAA10930). All of the active regions had concentrated, 
more or less balanced flux of both polarities and long tight (strong) PILs. 

4.4. Flux cancellations versus stepwise flux changes 

Do the flux cancellation events correlate with stepwise flux changes during the 
flares? Flux cancellation estimates during some flares show very large amount 
of flux loss (see Figure 5). We flnd features contributed to flux loss over the 
entire time interval aroimd a flare and only over during-fiare time. From them, 
we count number of features which involve stepwise fiux change. We do the 
comparison for the 45 flares where at least one feature with stepwise change 
was found. Table 13 shows that only small portion of features involved in flux 
cancellation showed stepwise change in their fluxes. Thus we find no relationship 
between these two types of event. The lack of relationship between the stepwise 
field changes and the cancellation suggests that the photospheric cancellation 
process is independent of the abrupt flare- related field changes. 



Table 13. Number of features with and without stepwise fiux 
changes contributed to flux cancellations over entire time interval 
around a flare and over duriiig-flare time. 



Number M-flares 


X-flares 




of features During flare Total 


During flare 


Total 


Stepwise flux change 2 


4 


13 


No step in flux 19 47 


37 


79 



5. Discussion 

We have analyzed magnetic flux changes during 77 major flares. We also com- 
puted flux cancellations associated with these flares. The results of this work can 
be summarized as follows: 
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We found that magnetic flux decreases occurring during the flares are signif- 
icantly more numerous than increases. The results are consistent with magnetic 
implosion and collapsing loop structure proposed by [Hudson, Fisher and Welsch 



( 2008 ) . Longitudinal field changes have larger values in the regions close to the 
limb, indicating that their field component is nearly horizontal. X-class fiares 
show larger values of the flux change and stronger bias towards flux decreases 
in comparison with M-class flares. 

Both unsigned and signed flux changes show overall increasing trend with 
flare strength. We find a linear correlation of the flux changes with GOES peak 
X-ray flux dominated by X-class flares in limb locations, in agreement with the 
study of jPetrie and Sudol| ( |2010 ) . 

We identified episodes of flux cancellation and estimated the amount of flux 
cancelled during three intervals around each flare; before flare onset, during flare 
and after flare end. We find that the amount of flux cancelled was greater over 
the time during the flare than before and after the flare, in general, with large 
contribution to the tendency from X-class flares. Some of the strongest X-class 
flares showed significantly greater fiux loss during the flare time. M-class flares 
did not show a significant bias in amount of flux cancelled during neither of the 
time frame. 

We note that the flux loss is approximately equally correlated with the number 
of canceling features and the total flux in these features that leads to the conclu- 
sion about approximately equal importance of the field strength and complexity 
in the flux cancellation processes during flaring activity. Increasing trend of the 
flux losses as a function of flare power is mostly explained by this correlation for 
the disk-center X-class flares, but not for the limb X-class flares. A less significant 
correlation of the flux losses with total flux in active region is also found for the 
disk-center locations, dominated by the X-class flares. 

The flux loss is correlated with the peak X-ray flux for both disk-center and 
limb X-class flares. Thus, correlation of both the flux losses and the flux changes 
with the flare strength is dominated by X-class fiares. 

We find that strong magnetic fields at PILs are associated with the occurrence 
of large flares in agreement with Schrijver (20071. We note that correlation with 
flare strength was higher for the flux loss than for the flux in canceling features. 

We did flnd contraction of the fleld towards neutral line in the form of flux 
cancellation events but we found no correlation between flux cancellation and the 
abrupt, stepwise photospheric fleld changes associated with coronal field implo- 
sion during flares. We speculate that the stepwise changes and fiux cancellation 
could be not related in general as the stepwise changes are caused by coronal 
magnetic restructuring and flux cancellation by photospheric dynamics. 

For future study, vector data is needed for a better interpretation of the longi- 
tudinal field changes results. High-cadence, high-sensitivity full vector data from 
SOLIS and HMI will be very useful. Also combined analysis of chromospheric and 
coronal data is needed to clarify the following: Are the canceling features part of 
the emerging U-loops or submerging i7-loops? Are the cancellations associated 
with reconnection? Do the reconnection occur above or below the photosphere? 
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